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ABSTRACT 

We present observations of far-infrared (50-200 pm) OH and H^O emission of the disk around the Herbig Ae star HD 163296 obtained 
with Herschel/PACS in the context of the DIGIT key program. In addition to strong [0 1] emission, a number of OH doublets and a 
few weak highly excited lines of H2O are detected. The presence of warm H2O in this Herbig disk is confirmed by a line stacking 
analysis, enabled by the full PACS spectral scan, and by lines seen in Spitzer data. The line fluxes are analyzed using an LTE slab 
model including line opacity. The H2O column density is 10 14 - 10 15 cirT 2 , and the excitation temperature is 200-300 K implying 
warm gas with a density n > 10 5 cm 4 . For OH we find N mot of 10 14 - 10 15 cnr 2 and T cx ~ 300-500 K. For both species we find an 
emitting region of r ~ 15 - 20 AU from the star. We argue that the molecular emission arises from the protoplanetary disk rather 
than from an outflow. This far-infrared detection of both H 2 and OH contrasts with near- and mid-infrared observations, which 
have generally found a lack of water in the inner disk around Herbig AeBe stars due to strong photodissociation of H2O. Given the 
similarity in column density and emitting region, OH and H 2 emission seems to arise from an upper layer of the disk atmosphere of 
HD 163296, probing a new reservoir of water. The slightly lower temperature of H 2 compared to OH suggests a vertical stratification 
of the molecular gas with OH located higher and H 2 deeper in the disk, consistent with thermo-chemical models. 
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1. Introduction 

Water is a key molecule for the chemical and physical evolution 
of protoplanetary disks. Together with O and OH, it forms the 
main reservoir of oxygen. The formation of water ice layers on 
dust grains may improve their sticking behaviour and thereby 
help the coagulation process towards larger particles that ulti- 
mately leads to the formation of planetesimals and planets. The 
growth of icy grains is also likely involved in the delivery of 
water to planets. Recent observations at near- (~ 3 pm) and mid- 
infrared (~ 10-30^/m) wavelengths have revealed the common 
presence of hot (T ~ 500-1 00 K) water and OH vapor in the 
1 atmosphere of T Tauri disks dSalvk et al.ll2b08t ICarr"& Naiita 
! I2008t lPontoppidanetal.ll2010t iMandell et all 1201 2|) . In con- 
trast, disks around the more massive Herbig AeBe stars do not 
show hot water vapor emission and appear to be depleted in 
water molecules dMandell et al.l 12008; Pontoppi dan et al.ll2010t 
iFedele et al.ll201 ll) . A likely explanation is that hot water is pho- 
todissociated by the stronger ultraviolet radiation emitted by the 
Herbig stars in the regi ons of the disks (< a few AU) probed at 
these wavelengths (e.g. lFedele et al.l201 ll) . Indeed, in the case of 
the young eruptive star EX Lupi the H2O emissio n is variable as 
a con sequence of the changing UV radiation field (Banzatti et al. 
120121) . 
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However, cooler water may survive further out or deeper into 
these disks, but that region can only be probed by longer wave- 
length data. Herschel offers the opportunity to search for these 
water lines with high sensitivity. Det ections of water i n disks 
with Herschel have been reported bvlHogerheiide et alJ (1201 ll) 
using HIFI an dRi viere-M arichalar et al. (201 1) using PACS, but 
these observations refer only to disks around T Tauri stars. 
In this letter we report the detection of OH far-infrared emission 
lines and the signal of warm H2O toward th e Herbig Ae star HD 
163296 (A1V) at a distance of d = 118 pc (Ivan Leeuwenll2007l) . 
The star is isolated with no evidence of a s tellar companion and 
is surrounded by a well- s tudied disk (e.g., [M annings & Sargent 
119971: iGradv et al J 120011: flsella et al.ll2007l) . A bipolar microiet 
and a series of Herbig-Haro knots are obser ved at optical and UV 
wavelengths perpendicular to the disk (e.g., Was sell et al.l2 006). 
The disk has recently been modeled by Till ing et alJ d2012l) . who 
also report upper limits on selected OH and H2O lines from 
PACS data obtained in the GASPS Herschel key program. 



2. Observations and data reduction 

HD 16329 6 was observed on A pril 03 2011 with the PACS 
instrument (Poglits ch et alJl2 010) onboard the Herschel Space 
Observatory dPilbratt et alJl2010» as part of the DIGIT key pro- 
gram {KPOTjievansJ , PI: N. Evans). The target was observed 
in SED mode covering the wavelength range 50-220 pm with 



1 



Fedele et al.: Warm H 2 and OH in the HD 163296 disk 



E 



-5 
6.1.0 



1 


■PpiT 


ft* 1 


1 1 ■■■ 





63.2 63.3 
X (tim) 



,1 5 

a 



E 





t 










71.5 
X (fim) 



E 



.O" 




d <£ 




l (um) 



84.0 84.2 84.4 84.6 84.8 85.0 
I (fim) 



Fig. 1. PACS spectrum (continuum-subtracted) of selected lines. The (blue) dashed line indicates the r.m.s. of the baseline multiplied 
by 3. The presented spectrum has been smoothed (smooth width = 2 bins) for clarity. The red line is a Gaussian fit to the detected 
lines. 



R ~ 1000 - 3000 (obsid: 1342217819, 1342217829). The ob- 
servations were carried out in chopping/nodding mode with a 
chopping throw of 6'. The total on-source integration time is 
6176 seconds for the B2A (51-73 fim) and short Rl (70-105 fim) 
modules and 8360 seconds for the B2B (70-105 fim) and long Rl 
(1 40-220 ^m) modules. The data have been reduced with HIPE 
8.0.2489 with standard calibration files from level to level 2. 
The two nod positions were reduced separately (oversampling 
factor = 3, up-sampling factor = 1 to ensure that the noise in 
each spectral point is independent) and averaged after a flat-field 
correction. 

The spectrum has been extracted from the central spaxel 
(9.4" square) to optimize the signal-to-noise (5 /AO ratio. Due 
to the large pont spread function of the telescope, some flux 
leaks into the other spaxels of the PACS array. To recover the 
absolute flux level, we apply a correction factor using the spec- 
trum extracted from the central 9 spaxels (3x3 extraction): this is 
performed by fitting a 3 r£, -order polynomial to two spectra (cen- 
tral spaxel and 3x3), the conversion factor is the ratio between 
the 2 fits. Finally, the spectrum is scaled so that the spectrum 
matches the PACS photometry (from Meeus et al. in prepara- 
tion) at 70jum and 160yum. 

The line flux (fii ne ) was measured by fitting a Gaussian 
function and the uncertainty (cr) is given by the product 
S TDf ^6A FWHAiH where STD F (W irT 2 fim' 1 ) is the stan- 
dard deviation of the (local) spectrum, 5 A is the wavelength spac- 
ing of the bins (fim) and FWHM is the full width half maximum 
of the line (fim). 



Table 1. [Oi] and OH line fluxes. 
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Notes. Column FJjL reports the line flux predicted by the best-fit model. 
<n) Blended with oTTO 6 25 -5 14 . 



Table 2. H 2 Q line fluxes. 
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3. Results 

We clearly detect the strong [O i] 63 fim line as well as 5 OH 
far-infrared features above 3<x (i.e. having Fn ne /cr > 3, TableQ}. 
Spectra of selected lines are shown in Fig. [T] The OH lines are 
readily recognized because of their doublet pattern; only intra- 
ladder transitions, which have the largest Einstein-A coefficients, 
are found. In the case of the 2 Tl[/2 7/2-5/2 doublet at 71 fim only 
one of the two lines is detected, although the non-detection of the 
second line is hardly significant within the noise. Asymmetric 
line intensities of A-dqublets are predicted at high temperature 
(Offer & van Dishoeck 1992), but because of the noise this dou- 
blet is not considered in our further analysis. 

Three lines of H2O are detected slightly above 3<r (Tabled. 
The H2O 818 - 7(>7 line at 63.32/mi is seen not only in our data 

1 This formula comes directly from the error propagation of the sum 
E,(F, ), where F; is the flux of the ;-th spectral bin. 



Notes. Flux integrated over 5 bins centered at the expected line posi- 
tion. 



but also in the GASPS spectrum shown bv lTilling etalj d2012l) 
although they do not claim the detection. The H2O 4 2 3 - 3 12 line 
at 78. 74//m is seen here w ith a flux of 1.8 (+ 0.4) 10~ 17 Wirr 2 
while Tilling et al.l (1201 2|) report only a 3<x upper limit of 1.5 
10~ 17 WirT 2 . lMeeus et al.l (1201 2h recently claim a detection of 
far-infrared H 2 Q emission towards this source based on new 
GASPS data. [Pontoppi dan et all (1201 Oh also provide tentative 
detections of H2O lines in the mid-infrared Spitzer wavelength 
range. Table [2] summarizes our fluxes and includes the fluxes 
measured at the position of some (undetected) key H2O lines 
that are used later in the analysis. 

The detected lines have upper level energies over a wide 
range with EJk ~ 120 - 900 K (OH) and EJk ~ 400 - 1300 
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K (H2O). Most of the lines are detected in the blue part of the 
spectrum. 
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Fig. 2. Stacking of 54 H2O lines. Water emission is clearly de- 
tected in the stacked spectrum. The integrated signal is 5h,o = 

7cr. 



3.1. Confirmation of H 2 by line stacking 

Since only 3 H2O lines are marginally detected above 3cr, we 
have used the availability of the full DIGIT PACS spectrum to 
confirm the presence of warm water in this disk through a stack- 
ing analysis. Line stacking is commonly done in extra-galactic 
surveys to detect th e faint emission lin es from the outer regions 
of galaxies (e.g.. ISchruba et all 120 1 ll) . Warm water has many 
lines spread throughout the far-infrared wavelength region which 
can be used for this purpose. In this work, we stacked spectra 
centered at the location of different H2O lines based on the far- 
infrared lines d etected with PACS t owards the protostar NGC 
1333 IRAS 4B dHerczeg et al.ll2012l) . The 95-100 fan range is 
excluded because of spectral leakage (produced by overlap of 
grating orders). Blended lines are excluded from this analysis 
and OH and [0 1] lines are masked. The remaining number of 
H2O lines available for the analysis is 54. The stacked spectrum 
is the weighted average of 54 spectra, each of which is 100 bins 

where 



wide centered at the position of a water line: F 

Fj is the (continuum-subtracted) spectrum centered at the y'-th 
water line and Wj is the weight of the line. The weight corre- 
sponds to STDj 1 , where STDj is the standard deviation of the 
continuum-subtracted spectrum Fj. The lines are stacked by bin 
because the spectral resolution in velocity space varies but is ap- 
proximately constant in bins. 

The stacked spectrum is shown in Fig. [2] The warm H2O 
signal is clearly detected and centered on the central bin. The 
integrated H2O signal is 7 times its uncertainty. The false alarm 
probability (FAP), i.e. the probability to detect a 7cr signal by 
stacking random portions of the PACS spectrum, is < 0.03% 
based on 10,000 randomized tests (see AppendixIBl. 



10 



10 17 k 



10 16 k 



10 15 lr 



10 14 ^ 



18 



10 



a 10 r 

£ io lf > 

c 

o 

a 10 15 k 

I 10 14 k 




100 



1000 



Excitation Temperature (K) 



Fig. 3. Contours of the reduced x 1 f° r the slab/LTE model for 
H2O (upper panel) and OH (lower panel). The value of the re- 
duced x 2 is overplotted. The (red) dashed lines provide the ra- 
dius of the emitting region (10, 20, 30 AU). The star indicates 
the location of the minimum^ 2 . 



This analysis confirms the presence of warm H2O in the 
PACS spectrum of HD 163296. Stacking H2O lines separately in 
spectra from the two nod positions also yields > 3<x detections. 
The H2O signal is only detected in the central spaxel and not in 
off-source spaxels. These last two tests exclude the contamina- 
tion from an extended and/or off-source emission and confirm 
that the H2O lines detected in the PACS spectrum are associated 
with HD 163296. 



4. Analysis 

In this section we analyze the OH and H2O excitation using a 
uniform slab of gas in local thermal equilibrium (LTE) including 
the effect of line opacity (see Appendix [A] for details). The lim- 
ited number of lines and their large uncertainties do not warrant 
a more sophisticated non-LTE treatment. The analysis is based 
on the data in Tables[T]and|2] i ncluding the mid-in frared lines de- 
tected with SpitzerflRS dPontoppidan et al] l2010). Free parame- 
ters of the model are the excitation temperature T ex (K) and the 
molecular column density N mo i (cm -2 ). The size of the emitting 
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Fig. 4. Spitzer/IRS spectrum (continuum-subtracted) of selected 
lines. The red and blue lines are the best model for H2O and OH 
respectively. 

region, given by its radius r, is not a free parameter since it can 
be determined uniquely for every given combination of T ex and 
A^moi. Comparison between models and data is done based on 
reduced^ 2 values. 

The range of models that yields an acceptable agreement is 
shown in Fig. [3] Overplotted are contours for the radius r. For 
H2O, the data are best fitted (lcr, p = 68.3 %) by models with T ex 
~ 200 - 350 K and N mo i ~ 10 14 - 10 16 cnr 2 and r ~ 15 - 20 AU. 
For OH, the data are best fitted for yV mo i ~ 10 14 - 10 15 cirT 2 , 
T ex ~ 300 - 500 K and r ~ 20 AU. The Spitzer/IRS spectrum of 
selected lines is shown in Fig.|4]along with the best model. 

Further constraints to the H2O column density and tempera- 
ture come from individual line flux ratios. In particular, the ratio 
between far- and mid-infrared lines (e.g. 7 07 - bib/lbi - 652) 
constrains the column density to > 10 14 cirT 2 . On the other 
hand, the ratio 8 is ~ 7m/8o$ — 7 \i constrain the upper limit of 
the H2O column density to 2 x 10 cm -2 . We also note that 
models with higher column density produce several H2O lines 
which are not detected in the PACS and IRS spectra. The ratio 
818 - 7o7/9o9 - 8 18 constrains the temperature to < 300 K. 

5. Discussion 

The primary result of this Letter is a detected signal of H2O to- 
ward a Herbig star, in addition to OH. The H2O emitting region 
is found to be 15-20 AU, demonstrating that H2O can survive 
the UV radiation further away from the star while it is likely 
photodissociated in the inner part of the disk. 

Given tha t a bipolar microjet is known to be associated with 
HD 163296 dWassell et alj 12006). the question arises whether 
the far-infrared molecular line emission presented here indeed 
arises from the disk or whether it comes from such a jet. There 
are several arguments in favor of the disk. First, we note that 
HD 163296 is isolated and not associated with a molecular 
cloud. No e vidence of a mo lecular outflow has been reported 
to date (e.g. iBae et alj|201l]) . Second, the spectrally resolved 
CO 7 = 3-2 line in the sub-millimeter shows the ch aracteris- 
tic double-peaked profi le for gas in Keplerian rotation dThi et alJ 
I200U iDent et alj|2005h . At much shorter wavelengths (4.7 pxa), 
the CO ro-vibrational em ission lines are al so characterized by 
a double-peaked profile (ISalvket al.ll201ll) . Thus, there is no 
hint of any significant small- or large-scale molecular outflow 
in these data that could dominate the PACS emission. Third, the 



PACS data show no evidence for extended/off-source emission 
beyond the central spaxel, not even for the strong [0 1] 63 pm 
line, which places the warm H2O within 500 AU of the central 
star. 

The inferred OH and H2O excitation temperatures of several 
hundred K indicate warm emitting regions. The high critical den- 
sities of the H2O lines, n c > 10 7 cm 4 , implies t hat the density 
of the gas should also be high (n > 10 5 cm~ 3 , e.g. lHerczeg et al.l 
12012b . These conditions and the arguments above suggest that 
the OH and H2O emission arises from the atmosphere of the disk 
associated with HD 163296 at radial distances >10 AU from the 
star. 

Assuming that the OH/H2O far-infrared lines are emitted by 
the disk, which zone does this emission trace? Models of the 
water chemistry in Herbig disks suggest at least three chem- 
ically distinct zones (e. g.. IWoitke et aU 120091 iGlassgold et al.l 
2009: IWalsh et al.l20 10. 2012: IVasvunin et al.ll201 ltlNaiita et all 
1201 ll) : (i) an inner disk water reservoir (<, few AU) with a chem- 
istry close to LTE; (ii) a cold water belt at large distances (>50 
AU) where gaseous H2O results primarily from photodesorption 
of water ice; and (iii) a hot water layer at intermediate distances 
of 1-30 AU and at medium heights with water formation driven 
by high temperature neutral-neutral reactions. The derived pa- 
rameters for our OH and H2O lines are consistent with zone 
(iii) (see also Tilling et al.ll2012l) : zone (i) is probed by the near- 
infrared data and zone (ii) can be targeted by HIFI observations 
of low-/ lines. Thus, the PACS data reveal a new water reservoir 
in disks. 



6. Conclusions 

We have presented new Herschel/PACS observations of the disk 
around the Herbig Ae star HD 163296. We obtain detections of 
far-infrared lines of warm OH and H2O toward a Herbig star. 
The presence of warm H2O is confirmed by a line stacking anal- 
ysis {lcr detection) enabled by the full PACS spectral scan. The 
LTE slab model analysis including optical depth effects indicates 
emission from the intermediate radii of the disk. Combined with 
near-infrared and sub-millimeter data, the oxygen chemistry can 
now be probed over the entire disk range. 
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Appendix A: Slab model 

For an optically thin line from a point-like source, the flux can 
be written by 



» —E u ]kT 

Fui = d£l s ■ Iui = d£l s —r-A ul N moi 8 " < L, rr . 

An Q(T) 



(A.l) 



with the solid angle of the source dQ. s , the line frequency 
v u i, the Einstein-A coefficient A„/, the molecular column density 
N mo i, the statistical weight of the upper level g u , the energy of the 
upper level E u and the partition function Q(T). The solid angle 
of the emitting region can be written as d£l s = nr 2 /d 2 , with the 
radius of the emitting region r and a distance of d — 1 18 pc to 
HD 163296. For an optically thick line, the integrated intensity 
is obtained from 



/„, = Av^B v „,(r ex )(l- e - T »'). 
c 

with the opacity at the line center 

A„/c 3 



Tul = 



Sttv^Av 



Ni— -N u 

8i 



(A.2) 



(A.3) 



The (thermal) width of the lines is assumed to be Av ~ 1 km 
s _1 , appropriate for gas at several hundred K and we assume a 
simple square like line p rofile as e.g. used in the RADEX code 
(Ivan der Tak et alj|2007l) . 
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Fig. B.l. Distribution of S jcr after 10,000 random stackings of 
54 parts of the PACS spectrum of HD 163296. The red line 
shows the Gaussian fit. The arrow indicates the location of the 
H 2 Q signal (Fig. |2|. 



Appendix B: False alarm probability of water 
detection in the stacked spectrum 

We performed a simulation to measure the probability to detect a 
signal with an integrated value, S > Icr. This provides the false 
alarm probability (FAP) of the detection in the stacked spectrum. 
We performed 10,000 random stackings of 54 (equal to the num- 
ber of water lines) parts of the PACS spectrum of HD 163296. 
After 10,000 iterations we measured the distribution of the ratio 
of the integrated signal over its uncertainty (measured as in Sec. 
13. 11 1. We mask the bins containing H2O, OH and [Oi] emission. 
Fig. IB. II shows the distribution of S jcr. The distribution is well 
fitted by a Gauss function (red line, x 1 = 0.03), centered (as ex- 
pected) at S l<x = (i.e. an equal number of positive and negative 
peaks). The number of occurrences with S jcr > 7 is < 3, which 
corresponds to FAP < 0.03% according to Bayesian statistics. 



